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I.   INTRODUCTION 

The  removal  of  heat  from  electronic  equipment  to  assure 
reliable  operation  has  been  a  problem  plaguing  the  electronic 
industry  for  many  years.   As  a  result,  many  researchers  have  - 
attempted  to  develop  specific  correlations  to  enable  designers  to 
attempt  to  predict  possible  solutions  to  this  problem. 

This  research  concerns  air  flowing  across  cylindrical 
elements  mounted  in  an  array  on  one  wall  of  a  vertical  adiabatic 
channel.   This  configuration  is  much  the  same  as  an  array  of 
transistors  in  the  T05  configuration  arranged  on  stacks  of 
parallel  circuit  boards  in  various  types  of  electronic  equipment. 

A.   THEORY  OF  VERTICAL  CHANNELS  FORMED  BY  CIRCUIT  BOARDS 

Many  electronic  systems  consist  partially  of  arrays  of 
printed  circuit  boards  supporting  a  fairly  dense  array  of 
electronic  elements  on  one  side  as  shown  in  Figure  1.1.   These 
circuit  boards  can  be  arranged  so  that  the  side  with  the  elements 
(front  side)  faces  the  relatively  smooth  backside  of  another. 
When  a  series  of  circuit  boards  are  arranged  in  such  a  fashion 
(see  Figure  1.2),  long  vertical  channels  in  which  cooling  air  can 
be  supplied  by  forced  convection  can  result.   During  operation, 
the  elements  such  as  flatpacks,  dips  resistors  and  transistors 
generate  heat,  the  majority  of  which  is  dissipated  to  the 
surroundings  with  little  transmitted  back  through  the  circuit 
board.   As  a  result,  it  appears  as  if  all  of  the  heat  is 
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Figure  1.1    Modern  Circuit  Board 
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Figure  1.2    Modular  Arrangement  of  Circuit  Board 
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transmitted  by  one  side  of  the  channel  or  one  side  of  each 
circuit  board  to  the  air  stream  or  air  layer  between  the  boards. 

Complexity  is  added  to  the  system  when  elements  serve 
different  purposes  and,  as  a  result,  possess  different  sizes, 
shapes  and  heat  dissipating  characteristics. 

The  process  of  heat  removal  from  electronic  systems  creates  a 
very  comlex  design  problem  governed  by  such  factors  as  circuit 
board/element  arrangement,  spacing,  heating  characteristics  and 
operating  environment.   Moreover,  the  ultimate  goal  of  any 
thermal  analysis  and  the  prediction  of  operating  temperatures  is 
the  assessment  of  reliability.   Reliability  is  the  probability 
that  an  electronic  component/system  will  operate  properly  without 
failure  for  its  intended  lifetime  under  its  environmental 
operating  conditions.   [Ref.  1] 

For  many  years,  there  have  been  numerous  attempts  to  derive 
techniques  to  aid  designers  in  this  field.   Unfortunately,  the 
methods  developed  are  often  too  complex,  inadequate  or  restricted 
to  specific  conditions. 

Early  work  with  natural  convection  between  parallel  arrays  of 
vertical  circuit  boards  was  performed  by  Elenbaas  [Ref.  2], 
Bodoia  and  Osterle  [Ref.  3]  and  Engel  and  Muller  [Ref.  4],   In 
1972,  Aung  et  al .  [Ref.  5]  studied  the  effects  of  circuit  board 
supports  and  elements  in  the  field  of  flow  for  wide  channel 
spacings.   Aung  went  on  to  investigate  uniform  heat  flux  on  two 
channel  walls  in  which  he  subsequently  studied  the  effects  of 
baffles  and  circuit  board  staggering.   [Ref.  6]   There  is 
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considerable  literature  describing  the  efforts  that  have  been 
made  to  attempt  to  model  circuit  channel  cooling  by  natural 
convection  and  forced  convection. 

For  example,  Kelleher  utilized  a  thermal  network  approach  to 
derive  a  more  representative  model  of  forced  convection  for  o.ne 
specific  geometry  with  equal  heat  releases  [Ref.  7].   Other 
experimental  or  semi -experimental  approaches  were  performed  by 
Tamura  [Ref.  8],  Finch  [Ref.  9],  Laermer  [Ref.  10],  Marto  [Ref. 
11]  and  Parkes  and  Preston  [Ref.  12].   Most  of  this  work  is 
restricted  to  specific  conditions. 

B.  THEORY  OF  SURFACE  ROUGHNESS  GEOMETRIES 

Surface  roughnesses  have  been  classified  by  several  different 
techniques  of  which  the  most  commonly  used  are  the  k  and  d  type 
roughness  classifications  consisting  of  2  and  3  dimension  classes 
which  are  further  divided  into  artificial  and  commercial 
roughness  subclasses.   Although  many  researchers  allude  to  these 
classifications,  only  a  limited  number  actually  use  this 
terminology  in  agreement.   It  is  further  noted  that  in  the  area 
of  heat- transfer,  the  k-  and  d-  type  roughnesses  seem  irrelevant 
because  the  same  heat  transfer  roughness  parameter  (g+) 
correlation  is  proposed  for  all  types  of  roughnesses. 

C.  THEORY  OF  HEAT  TRANSFER  FROM  ROUGH  SURFACES 

Earlier  researchers  in  this  field  seem  to  agree  that  the 
heat- transfer  roughness  parameter  (q+),  commonly  referred  to  as 
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the  inverse  of  the  sublayer  Stanton  number,  takes  the  functional 

form : 

P         q 
9+   =    1/Stsubiayer   =    C(e+)    (Pr) 

where   e+    is   the    roughness   Reynolds    number   based    on   an    arbitrary 
characteristic    height   of    the    roughness    (e).      The   different   va.lues 
of    p    and   q    fall    into    two   categories: 

1.  The   assumption   of    similarity    near   the   wall    which    results 
in   value    of    p    =    .2    and    q    =    .44    or    for    the    eddy-di f f usi vi ty 
models    p    =    .25    and    q    =    .66    [Refs.    13    and    14]. 

2.  The   assumption    of    a    breakdown    in   close    to    the   wall 
similarity    results    in    values    of    p    =    .45    and    q    =    .57    or 
p    =    .5    and    q    =    .66    for    the    eddy    diffusivity   models 
[Refs.    15   and    14]. 

With    regard    to    this    study,    the    foregoing    is    only    of    academic 

interest   because   the   geometry    here    is    defined    as    an    in-line    array 

of   cylindrical    elements. 
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II.   EQUIPMENT  DESIGN  AND  INSTRUMENTATION 

The  equipment  utilized  in  this  research  was  either 
manufactured  by  the  Engineering  Machine  Shop  or  supplied  by  the 
Naval  Engineering  Department  of  the  Naval  Postgraduate  School., 
Monterey,  California. 

A.   WIND  TUNNEL 

The  wind  tunnel  (Figure  2.1)  is  constructed  of  1/8" 
plexiglass  with  the  exception  of  two  1/4"  access  plates  and 
various  reinforcement  members  mounted  on  the  exterior  so  as  to 
ensure  smooth  air  flow  on  the  tunnel  interior. 

The  test  section  of  the  tunnel  has  a  12"  by  4"  cross  section 
and  a  height  of  18"  with  access  covers  on  both  faces  to 
facilitate  instrumentation  and  support  for  the  test  surfaces. 

The  tunnel  intake  is  belled  with  an  18"  radius  to  produce  an 
orifice  that  has  a  40"  by  12"  cross  section.   An  adjustable  leg 
is  attached  at  each  corner  to  serve  as  a  support  for  the  tunnel, 
to  level  the  tunnel  and  to  provide  an  8"  clearance  off  the  floor. 

The  upper  end  of  the  tunnel  is  reduced  to  a  section  of  2"  PVC 
ducting  which  is  connected  to  a  550  cfm  turbo  fan.   Air  flow 
through  the  tunnel  is  controlled  by  a  series  of  valves  connected 
in  line  with  the  fan. 

A  1/8"  plexiglass  splitter  plate  was  mounted  1/2"  off  of  the 
access  plate  facing  the  test  surface  in  order  to  reduce  the 
effects  of  the  boundary  layer  formed  in  the  tunnel.   This 
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Figure  2.1    Photograph  of  Wind  Tunnel 
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distance  was  pre-de termi ned  through  the  use  of  computer  software 
[Ref.  16]  and  the  results  are    displayed  in  Figure  2.2.   A  wooden 
frame  box  (4'x4'x  81)  was  constructed  around  the  wind  tunnel  to 
elliminate  the  effects  of  incident  radiation  on  the  cylinders 
from  the  laboratory  environment  and  the  box  also  served  to  reduce 
any  transient  air  currents  created  by  movement  in  the  room  (i.e., 
doors  opening  and  closing). 

B  .   TEST  SURFACE 

The  design  of  the  test  surface  used  in  this  research  was 
derived  from  the  previous  work  of  Arvizu  and  Moffat  [Ref.  17]  and 
Ortega  [Ref.  18].  It  was  constructed  of  a  12"  x  12.5"  x  1/8" 
sheet  of  plexiglass  (for  support)  with  a  12"  x  12.5"  x  3/8"  sheet 
of  balsa  wood  (thermal  boundary)  cemented  to  it.   The  leading 
edge  of  this  (balsa  wood-plexiglass)  splitter  plate  was  cut  at  a 
35°  angle  to  form  a  sharp  leading  edge  to  reduce  the  effects  of 
leading  edge  separation.   The  surface  (balsa  wood)  of  the 
splitter  plate  was  then  sanded  and  painted  three  times  with  white 
enamel  paint  and  sanded  a  fourth  time  to  ensure  a  smooth  test 
surface.   Small  (1/16"  diameter)  holes  were  drilled  through  the 
plate  to  permit  the  installation  of  resistor  leads  and 
thermocouples  into  the  cylindrical  elements  mounted  on  the  balsa 
wood.   Three  mounting  screws  were  attached  to  the  plexiglass  side 
of  the  splitter  plate  and  to  the  access  plate  to  facilitate  the 
positioning  of  the  test  surface.   A  tail  fin  mounted  to  the  tail 
end  of  the  splitter  plate  was  used  to  regulate  the  air  flow  on 
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both  side  of  the  splitter  plate.   A  fourth  adjusting  screw  was 

attached  to  this  fin  and  through  the  access  plate 

(Figure  2.3). 

C.   CYLINDRICAL  ELEMENTS 

The  cylindrical  elements  used  in  this  research  were 
manufactured  from  3/4"  round  aluminum  stock  and  cut  to  a  height 
of  1.0  cm.   A  7/16"  diameter  hole  was  cut  to  a  depth  of  1/4"  in 
the  center  of  the  element  and  two  additional  holes  (1/8" 
diameter,  1/4  deep)  were  cut  opposite  each  other  to  permit  the 
installation  of  the  56  ohm,  1/2  watt  resistors.   Another  hole 
(1/32"  diameter,  3/8"  deep)  was  drilled  to  position  the 
thermocouple  (Figure  2.4). 

The  resistors  were  epoxied  into  forty  two  of  the  cylinders 
and  six  more  were  cut  in  half  so  that  an  array  of  7x6  cylinders 
could  be  mounted  evenly  on  the  splitter  plate  and  the  12  half 
cylinders  could  be  arranged  along  the  edges  of  the  plates  (Figure 
2.5)  . 

The  resistor  (56  ohm,  1/2  watt),  prior  to  insertion  into  the 
cylinders,  had  the  leads  cut  and  smaller  leads  (2841/1,  #30  solid 
spc.)  attached  so  as  to  reduce  heat  losses  by  conduction  through 
them. 

The  cylinders  were  mounted  on  the  test  surface  in  an  in-line 
array  with  a  separation  of  one  cylinder  diameter  between 
cylinders.   This  was  done  both  vertically  and  horizontally. 
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Figure  2.3   Sketch  of  Test  Surface  Configuration. 
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Top  Row:   Assorted  Views  of  Cylinders 
Bottom  Row:   T05  Transistor  and  a  nickel 


Figure  2.4   Photograph  of  Cylinders. 
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Figure  2.5    Photograph  of  Cylinders  Arranged  in  an  Array 
on  the  Test  Surface. 
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D.   INSTRUMENTATION 

1 .  Temperature  Measurements 

Type  T,  copper-cons tantan  thermocouples  were  utilized  for 
the  determination  of  the  temperature  of  the  heated  elements  and 
the  inlet  and  outlet  temperatures  of  the  wind  tunnel  near  the 
test  surface.   The  thermocouples  were  connected  to  seven  scanner 
boards  in  an  autodata  nine  and  calibrated  to  an  accuracy  of 
+_0.2°F  (Appendix  A.)   The  thermocouples  were  then  routed  through 
an  access  plate  in  the  tunnel  through  the  splitter  plate  and  into 
each  of  the  elements. 

2 .  Velocity  Measurements 

The  velocity  measurements  were  obtained  by  the  use  of  a 
micro  manometer  and  a  pitot-static  tube.   Changes  in  the  dynamic 
and  static  pressure  were  measured  by  the  micro  manometer  and 
using  Barnoulli's  equation,  velocity  was  determined.   The  pitot- 
static  tube  was  placed  between  the  splitter  plate  and  the  test 
surface  splitter  plate  and  used  to  determine  the  approach 
velocity  between  them. 

A  TSI  (Model  1050/1050AA)  hot  wire  anemometer  was  used  to 
determine  the  velocity  profile  prior  to  entry  into  the  test  area. 
The  hot  wire  anemometer  was  calibrated  with  the  micro  manometer 
and  used  only  to  ensure  that  the  approach  velocity  was  uniform. 

Accuracy  of  the  mirco  manometer  was  determined  to  be 
+_.0005  inches  of  water  and  resulted  in  an  error  of  +0.2  ft/sec  or 
+7cm/s.   Because  the  micro  manometer  was  utilized  for  the 
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calibration  of  the  hot  wire  anemometer,  more  accuracy  might 
result  if  the  micro  manometer  was  used  to  determine  the  velocity 
between  the  splitter  plates. 

3 .  Power  Measurements  and  Control 

A  LAMBDA  (60  volt,  10  amp)  Regulated  Power  Supply  (Model 
LK  345AFM)  was  used  to  power  to  the  56  ohm,  1/2  watt  resistors 
mounted  in  the  forty  two  cylindrical  elements. 

A  switch  board  was  set  up  to  allow  power  to  be  delivered 
to  individual  elements  in  either  of  two  columns,  the  third  or  the 
fourth  columns.   It  was  also  set  up  in  such  a  manner  that  each  or 
all  rows  could  be  energized. 

A  Hewlett  Packard  digital  multimeter  (3466A)  was  used  to 
measure  the  voltage  drop  across  the  resistors  and  a  Fluke  8060A 
multimeter  was  used  to  measure  the  current  through  the  circuit. 
The  power  delivered  is  the  product  of  current  and  voltage. 

4 .  Data  Acquisition 

The  Autodata  Nine  mentioned  above  is  a  self  contained 
recordinq  device  with  up  to  100  channels  available  for  both 
temperature  and  voltage  measurements.   During  this  research,  71 
channels  were  utilized,  one  for  voltage  measurement  and  the  rest 
for  temperature  measurements  in  degrees  Fahrenheit.   (Figure  2.6) 

E.   EXPERIMENTAL  PROCEDURE 

For  this  research,  data  was  recorded  at  three  different 
channel  height  to  cylinder  height  ratios  (1.0,  2.3  and  4.6)  for 
two  different  approach  velocities  (1.5  m/s  and  3.0  m/s).   The 


26 


Figure  2.6    Photograph  of  Autodata  Nine 
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data  recorded  consisted  of  the  voltage  supplied  to  the  circuit, 
the  current  through  the  circuit  and  the  temperatures  of  the  test 
elements  as  well  as  the  inlet  and  exhaust  temperatures  of  the 
airstream.   For  each  channel  height  and  approach  velocity,  data 
was  gathered  after  a  steady  state  condition  was  reached,  for 
each  element  of  a  column  energized  one  at  a  time  and  for  the 
first  four  rows  and  for  all  six  rows  energized  all  at  once. 
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III.       DATA    REDUCTION 

A.  DATA  TAKEN 

Inlet  Free  Stream  Temperature,  T   (°F) 
Cylinder  Surface  temperature,  Ts,0  (°F) 
Atmospheric  Pressure,  Patm  (inches  Hg.) 
Pitot-Static  Pressure  Difference,  AP  (inches  H2O) 
Power  Supply  Voltage,  V  (volts) 
Current  Drawn,  I  (amps) 
Ts>i  (°F) 

B.  SOUGHT  AFTER  PARAMETERS 

1 .   Dimension! ess  Temperature  Ratio 

The  Dimensionl ess  Tempera ture  Ratio,  0,  is  obtained  as  a 
function  of  cylinder  location,  channel  configuration  and 
freestream  velocity.   Here  e   is  defined  as 


A  Ts,i  "  T 


Ts,l  -  T 


where  Ts,i  is  the  surface  temperature  obtained  by  the  first 
passive  cylinder  directly  after  the  cylinder  heated  by  a 
resistor.   Ts,j  is  the  surface  temperature  obtained  by  individual 
passive  cylinders  after  and  in  the  same  column  as  the  cylinder 
heated  by  the  resistor  where  i  is  the  row  after  the  heated 
cylinder.   Note  that 

Ts  1  "  T~ 
91  =  y^^ j —  =  1.0  by  definition 

1     's,l  "  '  °° 
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2 .   Heat- transfer  Coefficient 

The  Heat-transfer  Coefficient,  h(W/m2  c),  is  found  and 
is  a  function  of  the  cylinder  location,  channel  configuration  and 
freestream  air  velocity. 

The  heat  transfer  coefficient  was  determined  as  follows: 


q      (W)  -  q       (W) 

h  =    input   losses 

A  (Ts,o  -  Too  ) 


Ts,0  and  T^  were  recorded  by  the  autodata  nine  and 
corrected  by  the  calibration  data  in  Appendix  A. 

A,  is  the  surface  area  of  the  cylinder  is  calculated  for 
channel/ cylinder  height  ratios  of  1.0,  2.3  and  4.6.   If  the 
channel /cyl i nder  height  ratio  is  equal  to  1.0,  A  =  udL  and  if  the 
channel /cyl  i  nder  height  ratio  is  greater  than  1.0,  A  =  tt  [  dL  + 
d2/4]. 

The  value  of  Qiosses  in  wa^ts  is  derived  from  a  parabolic 
temperature  dependent  approximation  utilized  by  Ortega  [Ref.  18] 
in  his  researach  with  aluminum  cubes.   This  approximation  is 
assumed  to  be  the  same  for  cylinders  and  is  used  here. 

Qiosses  (W)  =  "0,62927  x  10"3  +  0.5059  x  1 0"2 ( Ts , q-T^ ) 

-  0.023695  x  10-5  (TS,0-TJ2 

Ortega  also  determined  the  radiation  losses  to  be  an  order  of 
magnitude  smaller  than  those  given  above,  and  radiation  losses 
were  therefore  disregarded. 
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Qinput(w)  1S  tne  Product  of  power  supply  voltage  and 
current  drawn. 

Qinput  =  V(volts)  x  I(amps) 

3 .   Reynolds  Number 

The  Reynolds  number  is  developed  as  a  function  of  channel 
configuration  and  free  stream  air  velocity.   In  this  research, 
the  Reynolds  number  was  calculated  using  the  freestream  air 
velocity,  channel  height  and  viscosity  as  follows: 


Re   =  channel  height  x  velocity 
H  viscosity 


The  value  of  the  viscosity  is  obtained  by  a  linear 
interpolation  of  data  from  tables  of  viscosity  and  air 
temperature. 

4 .   Nusselt  Number 

The  Nusselt  number  is  prescribed  as  a  function  of  channel 
configuration.   Here,  the  Nusselt  number  was  calculated  using 
channel  height,  heat- transfer  coefficient  and  the  thermal 
conductivity  of  the  free  stream  air  as  follows: 


Nu   =  channel  height  x  heat-transfer  coefficient 

thermal  c  o  n  d  u  c  t  i  v  i  ty 


The  value  of  the  thermal  conductivity  was  obtained  by  a 
linear  interpolation  of  data  from  tables  of  thermal  conductivity 
and  air  temperature. 
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C.   TYPICAL  SET  OF  DATA 

Table  I  presents  a  typical  set  of  data  obtained.   These  data 
will  be  used  here  to  illustrate  the  data  reduction  procedure. 

TABLE  I 
DATA  SET 


patm  =  29.736  inches  Hg. 

AP  =  0.0051  inches  H20 

Channel/Cylinder  Height  Ratio  =  4.6 

Cylinder  Heated  -  row  1,  column  3 

Cylinder  Surface  Temperature,  Ts,q  =  206. 1°F 

Voltage  Supplied  =   14.380  volts. 

Current  Drawn     =   0.225  amps. 


Ts-1 

Ts>2 

Ts»3 

Ts>4 

Ts»5 
T 


80.3°F 
74.0°F 
72.2°F 
71.5°F 
70.8°F 
68.6°F 


D.   CALCULATIONS 

1 .   Dimension! ess  Temperature  Ratio, 

The  temperatures  for  the  cylinder  surface  were  recorded 
by  the  autodata  nine  and  corrected  by  use  of  the  thermocouple 
calibration  data  presented  in  Appendix  A. 

Ts,-j  (corrected)  =  m  x  Ts,-j  (recorded)  -  b 
m   Recorded    Jd    ts,  i  (Corrected) 

79.6  °F 


TS>1 


1.0026   (80.3)  -  0.9038 
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Ts,2  =  0.99994  (74.0) 
Ts,3  =  1.0018  (72.2) 
Ts,4   =   1,0018   (68.6) 

a  T  s  ,  i  -   f  a, 

6i   =    

TS»1   "   Too 


0.19287 
0.05645 
0.5980 


73.8      F 


72.3      F 


68.1 


1 

=    79.6    - 
79.6   - 

68.1      « 
68.1 

«      11.5      = 
11.5 

1.0    (b 

2 

=    73.8    - 
79.6   - 

68.1      = 
68.1 

5.7      = 
11.5 

0.4956 

3 

=    72.3    - 
79.6   - 

68.1      = 
68.1 

4.2      = 
11.5 

0.3652 

4 

=    71.2    - 
79.6   - 

68.1      = 
68.1 

3.1      = 
11.5 

0.2696 

5 

=    70.3    - 
79.6   - 

68.1      = 
68.1 

2.2      = 
11.5 

0.1913 

2.      Heat-Tansfer  Coefficient 


h    =    Q input   -    Qlosses 
A    (Ts,q    -    Too    ) 

Ts,0    =    1.0025    (206. 1°F)    -    0.4745    =    206. 1°F 

T^        =    1.0018    (68.6°F)       -    0.5980    =         68.1°F 

Ts,0    -    Too  =    206.1    -    68.1    =    138. 0°F    =    76.7°C 


(corrected) 
(corrected) 


Channel /cyl i nder   height    ratio   =    4.6 

therefore   A    =  tt  [(  0.019m)    (0.01m)    +    (  0  .019m)  2/4] 
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A    =    0.00088m2 

Qlosses    =    0.6927    x    10~3    +    0.5059    x    10"2  (  Ts  .q-T,,  ) 
-    0.23695    x    10-6(Ts>o-Too)2 

o 

where  Ts>0   -Too  is    in      C    and   Qlosses    1  s    ' n   watts. 

o 

f  s  » o"^°°  =    7  6.7    C 


Qlo 


sses 


=    0.373   watts 


Qinput   =    VxI    =    14.380    volts    x   0.225    amps    =    3.27   watts 

h    =    3.27    -    0.373        =    42.92    W/m2    °C 
0.00088    (76.7) 


3 .      Reynolds    Number 


ReH    =    HUro/  v 

H    =    4.6    (  .01m)    =    0.046m 

U^  =    1.5m/s 

=    38.02    x    10~6    -    109.8    x    10"9    (27.0 


-    T^   ) 


where   T         i  s    1 n      C 


T      =    20.0    C 

00 

=    38.02    x    10-6m2/s 

Re^    =    .046m    x    1.5m/s 

38.02    x    10"6m2/s 


=  1814.7 


4 .  Nusselt  Number 

NuH  =  Hh/k 

H  =  .046m        h  =  42.92  W/m2°C 

k  =  0.0407  -  76.8  x  10"6  (27.0-T  ) 
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where   T  ^     is    in      C 


k    =    40.162    x    10~3w/m°C 

NuH    =    .046m    x    42.92    W/m2°C 
40.162    x    10-3    W/m°C 


=    49.16 
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IV.   RESULTS 

The  experimental  measurements  that  were  obtained  lead  to  two 
different  subject  categories:   thermal  wake  effects  and 
heat- transfer  coefficient.   Table  II  presents  the  range  of  the 
parameters  covered  and  a  listing  of  the  data  is  tabulated  in 
Appendix  C. 

TABLE  II 
PARAMETERS  TESTED 

Thermal  Wake 


Channel  Height/Cylinder  Height  Ratio     1.0 
Downstream  Row  Number                0,1,2,5 
Vel  ocity  (m/s) 1.5,  3.0 


4.6 
0,1,2,5 
1.5,  3.0 


Heat- Transfer  Coefficient 


Channel  Height/Cylinder  Height  Ratio   1.0  2.3  4.6 

Heated  Cylinder  Row  Number          1,2,3,6  1,2,3,6  1,2,3,6 

Number  of  Rows  Heated                 4,6  4,6  4,6 

Velocity  (m/s) 1.5,3.0  1.5,3.0  1.5,3.0 

A.   THERMAL  WAKE 

The  thermal  wake  created  by  a  single  heated  cylinder  in  an 
array  is  influenced  by  the  channel  height  to  cylinder  height 
ratio,  the  mass  flow  rate  and  the  position  of  the  cylinder  in 
the  array. 

For  a  channel  height/cylinder  height  ratio  of  1.0,  there,  is 
no  free  stream,  and  the  flow  through  the  array  undergoes  periodic 
accelerations  and  decelerations.   It  therefore  seems  unlikely 
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that  a  channel  height/cylinder  height  ratio  of  4#6  woui(j  possess 
the  same  convective  means  for  determination  of  the  downstream 
wake.   However,  a  plot  of  6  as  a  function  of  downstream  row 
number  (Figure  4.1  and  4.2)  demonstrates  a  1/N  dependence  for 
both  velocities  at  both  channel  heights. 

Figure  4.3  shows  an  even  closer  resemblance  to  the  1/N  curve 
for  the  thermal  wake  created  by  a  heated  cylinder  located  in  a 
row  internal  to  the  matrix  of  cylinders. 

The  following  is  provided  for  the  understanding  of  the 
graphical  curves  found  in  the  legends  of  various  graphs  in 
Chapter  IV: 

Hx  .  xRxVx . x 

Hx.x   Channel  Height/Cylinder  Height  Ratio  x.x 
Rx     Row  that  the  heated  cylinder  is  in  x 
Vx.x    Approach  Velocity  x.x 

For  example,  H.4.6R1V3.0  denotes  data  gathered  at  a  channel 
height/cylinder  height  ratio  of  4.6,  where  the  cylinder  heated 
is  in  the  first  row  and  the  approach  velocity  is  3.0  m/s. 

B.   HEAT-TRANSFER  COEFFICIENT 

In  Figure  4.4  the  heat- transfer  coefficient,  h,  for  a  single 
heated  cylinder  in  the  first  row  is  plotted  against  channel 
height/cylinder  height  ratios  for  approach  velocities  of  1.5  m/s 
and  3.0  m/s.   It  can  be  seen  that  as  the  channel  height 
increases,  the  value  of  h  drops  off  fairly  rapidly  to  an 
asymototic  limit  at  a  channel  height/cylinder  height  ratio  of 
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4.6.   Air  flows  over  the  cylinders  instead  of  through  them  (the 
path  of  least  resistance). 

Figure  4.5  show  a  plot  of  h  against  the  channel  height  to 
cylinder  height  ratio  for  an  internal  cylinder  where  h  drops  off 
a  little  more  rapidly. 

Figure  4.6  plots  h  against  array  position  for  a  channel 
height/cylinder  height  ratio  of  1.0.   Here,  h  increases  with 
position  upstream.   This  correlates  to  the  increased  turbulence 
with  position  downstream. 

Figure  4.7  is  a  plot  of  h  against  array  position  for  a 
channel  height  to  cylinder  height  ratio  of  4.6.   For  this 
configuration  h  appears  to  increase  slightly  and  then  drop  off 
slightly  so  that  the  value  of  h  at  the  sixth  row  is  less  then 
that  of  the  first  row. 

Figure  4.8  plots  Nusselt  number  against  Reynolds  number  for 
all  three  channel  height  to  cylinder  height  ratios.   This 
composition  leads  to  a  correlation  of  the  form  Nu  =  C  Re^*5^ 
where  C  is  different  for  different  array  densities  as  noted  by 
Arvizu  and  Moffat  [Ref.  17]. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

In  this  research,  for  cylindrical  elements  mounted  in  a 
vertical  channel,  the  characteristic  curves  produced  closely 
resemble  those  for  a  similar  configuration  of  cubical  elements 
which  were  investigated  by  of  Arvizu  and  Moffat  [Ref.  17]. 

Although  there  is  room  for  more  data  to  assist  in  making  firm 
conclusions  about  the  flow  characteristics  of  air  flow  around 
cy 1 i  nders ,  the 

The  following  has  been  observed: 

1.  The  1/N  characteristic  curve  produced  by  the  thermal  wake 
of  a  heated  cylinder  on  an  adiabatic  wall  in  a  channel, 
may  allow  for  fairly  accurate  prediction  of  temperature 
increase  in  downstream  elements  in  an  array  for  various  air 
fl ow  rates. 

2.  Channel  spacings  and  airflow  rates  have  very  little  affect 
on  the  shape  of  the  characteristic  curve  produced  by  the 
thermal  wake  of  a  heated  cylinder  in  an  array. 

3.  The  position  of  a  heated  cylinder  in  an  array,  on  an 
adiabatic  wall  in  a  channel,  has  a  definite  affect  on  its 
heat  transfer  characteristics.   Elements  requiring  more 
cooling  should  be  placed  closer  to  the  freestream  air  inlet, 
This  would  reduce  effects  created  by  the  thermal  wake  of 
other  cyli nders. 

4.  Heat  transfer  in  an  array  of  cylindrical  elements  is 
mainly  influenced  by  the  air  flow  between  the  cylinders 
and  slightly  less  by  the  turbulence  in  the  air  flow. 

5.  As  the  channel  height  to  cylinder  height  increases,  the 
amount  of  heat  transfer  that  occurs  decreases  to  an 
asymototic  limit  near  a  channel /cyl i nder  height  ratio  of 
4.6. 

For  continued  work  in  this  field,  it  is  recommended  that,  a 

longer  test  surface  with  more  rows  of  elements  be  utilized  along 
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with  higher  velocities  and  various  geometric  shapes.   All  of  this 
can  be  accomplished  with  only  a  few  modifications  to  the 
equipment  at  hand. 
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APPENDIX    A. 
THERMOCOUPLE    CALIBRATION 
The   thermocouples   were    all    connected    to    a   wooden    framed 
switchboard   which   was   connected    by    10'    wire   leads    to    seven 
scanner   boards   mounted    in   the    autodata    nine   data    acquisition 
device.      This   combination   was    then   calibrated    as    a   complete   unit 
in   the   calibration   laboratory    against    five    known    temperatures: 

66.02°F,    106.14°F,    127.615°F,    151.20°F    and    200.42°F 
Thermocouple   readings    follow: 


1        1 

(6.3 

106.4 

127.9 

151.3 

200.5 

2            ( 

it. 3 

106.5 

128.0 

151.3 

200.6 

3            { 

>6.2 

106.5 

128.0 

151.2 

200.4 

« 

>6.J 

106.5 

128.0 

151.3 

200.6 

5 

,6.1 

106.4 

127.9 

151.2 

200.4 

C 

>6.2 

106.4 

127.9 

151.2 

200.4 

7 

.6.1 

106.2 

127.8 

151.1 

200.3 

8 

>6.2 

106.4 

127.8 

151.3 

200.4 

9 

,6.1 

106.2 

127.7 

151.0 

200.2 

10             i 

,6.2 

106.4 

127.8 

151.2 

200.4 

11             1 

,6.3 

106.9 

128.5 

151.8 

201.0 

12             1 

,6.8 

107.0 

128.5 

151.8 

200.9 

15             { 

>6.7 

106.8 

128.3 

151.6 

200.3 

14             1 

,6.7 

107.0 

128.4 

151.8 

201.0 

15             1 

>6.6 

106.8 

128.3 

151.6 

200.7 

16             1 

,6.7 

106.8 

128.3 

151.7 

200.8 

17             1 

,6.5 

106.7 

128.0 

151.4 

200.7 

18             1 

,6.7 

106.8 

128.2 

151.6 

200.8 

19             1 

(4.5 

106.6 

128.0 

151.4 

200.7 

20             ( 

,4.2 

106.8 

122.6 

152.4 

199.3 

21             1 

>6.2 

106.4 

128.0 

151.3 

200.4 

22             1 

.6.3 

106.5 

128.0 

151.3 

200.5 

23             1 

•  6.1 

106.3 

127.9 

151.2 

200.4 

24             ( 

,6.2 

106.5 

128.0 

151.3 

200.5 

25             1 

,6.0 

105.3 

127.8 

151.2 

200.3 

26             1 

,6.2 

106.5 

127.8 

151.2 

200.4 

27             1 

,6.0 

106.2 

127.1 

151.0 

200.2 

28             1 

,6.1 

106.2 

127.7 

151.1 

200.2 

29             1 

,6.0 

106.2 

127.7 

151.1 

200.3 

30             1 

,6.3 

106.5 

127.8 

151.2 

200.1 

31             1 

,6.3 

106.4 

127.8 

151.2 

200.4 

32             1 

,6.3 

106.5 

128.0 

151.4 

200.6 

33             1 

,6.2 

106.3 

127.9 

151.4 

200.6 

34             ( 

,6.3 

106.5 

128.0 

151.4 

200.6 

35 

66.3 

106.5 

127.9 

151.2 

200.4 

36 

66.3 

106.5 

128.0 

151.4 

200.5 

37 

66.3 

106.5 

127.9 

151.3 

200.5 

38 

66.4 

106.5 

128.0 

151.3 

200.6 

39 

66.3 

106.4 

127.9 

151.2 

200.4 

40 

66.6 

106.7 

128.2 

151.6 

200.7 

41 

66.5 

106.6 

128.1 

151.5 

200.7 

42 

66.6 

106.8 

128.4 

151.6 

200.9 

43 

66.6 

106.7 

128.2 

151.5 

200.8 

44 

66.7 

106.8 

128.4 

151.6 

200.8 

45 

66.5 

105.8 

128.2 

151.5 

200.7 

46 

66.6 

106.8 

128.3 

151.6 

200.8 

47 

66.5 

106.7 

128.2 

151.6 

200.7 

48 

66.7 

106.7 

128.2 

151.6 

200.7 

49 

66.6 

106.6 

128.0 

151.5 

200.7 

50 

66.4 

106.5 

128.1 

151.4 

200.5 

51 

66.3 

106.4 

127.9 

151.3 

200.6 

52 

66.4 

106.6 

128.2 

151.4 

200.7 

53 

66.3 

106.6 

128.1 

151.4 

200.6 

54 

66.4 

106.6 

128.2 

151.4 

200.7 

55 

66.4 

106.6 

128.0 

151.4 

200.6 

56 

66.4 

106.6 

128.3 

151.5 

200.7 

57 

66.4 

106.5 

128.0 

151.4 

200.6 

58 

66.6 

106.6 

128.2 

151.4 

200.6 

59 

66.5 

106.6 

128.1 

151.4 

200.6 

60 

66.2 

106.8 

128.0 

151.3 

200.4 

61 

66.2 

106.7 

128.0 

151.4 

200.5 

62 

66.3 

107.1 

128.2 

151.4 

200.6 

63 

66.3 

106.8 

128.0 

151.3 

200.6 

64 

66.3 

106.9 

128.2 

151.4 

200.6 

65 

66.3 

106.8 

128.1 

151.4 

200.6 

66 

66.3 

106.3 

128.1 

151.4 

200.6 

67 

66.4 

106.9 

128.3 

151.5 

200.7 

68 

66.3 

106.8 

128.1 

151.5 

200.7 
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A   least   squares    fit   program  was   developed    to   produce 
thermocouple   corrections   with   the   following    result: 
Tactual    =      M   x   Tmeasured   +   B 


THERMOCOUPLE 

1 

Ma 

THERMOCOUPLE 

2 

Ma 

THERMOCOUPLE 

3 

M» 

THERMOCOUPLE 

4 

rt« 

THERMOCOUPLE 

S 

Ma 

THERMOCOUPLE 

f 

M« 

THERMOCOUPLE 

7 

Ma 

THERMOCOUPLE 

8 

Ma 

THERMOCOUPLE 

9 

Ma 

THERMOCOUPLE 

10 

Ma 

THERMOCOUPLE 

11 

Ma 

THERMOCOUPLE 

12 

Ma 

THERMOCOUPLE 

IS 

Ma 

THERMOCOUPLE 

14 

Ma 

THERMOCOUPLE 

15 

Ma 

THERMOCOUPLE 

14 

Ma 

THERMOCOUPLE 

17 

Ma 

THERMOCOUPLE 

18 

Ma 

THERMOCOUPLE 

19 

Ma 

THERMOCOUPLE 

20 

Ma 

THERMOCOUPLE 

21 

Ma 

THERMOCOUPLE 

22 

Ma 

THERMOCOUPLE 

23 

Ma 

THERMOCOUPLE 

24 

Ma 

THERMOCOUPLE 

25 

Ma 

THERMOCOUPLE 

24 

Ma 

THERMOCOUPLE 

27 

Ma 

THERMOCOUPLE 

28 

Ma 

THERMOCOUPLE 

29 

Ma 

THERMOCOUPLE 

30 

Ma 

THERMOCOUPLE 

31 

Ma 

THERMOCOUPLE 

32 

Ma 

THERMOCOUPLE 

33 

Ma 

THERMOCOUPLE 

3« 

Ma 

1.0017233  Ba 

1.0012808  Ba 

1.00224*7  Ba 

1.0012808  Ba 

1.0013437  Ba 

1.001*844  Ba 

1.0017433  Ba 

1.0017471  Ba 

1.0024321  Ba 

1.001*540  Ba 

1.0017424  Ba 

1.0024770  B> 

1.0024207  Ba 

1.0013180  Ba 

1.0024844  Ba 

1.0024144  Ba 

1.0021400  Ba 

1.002595*  Ba 

1.0018*5*  B* 

0.9914704  Ba 

1.0018024  Ba 

1.0019808  Ba 

1.0011015  Ba 

1.0013380  Ba 

0.9984798  Ba 

1.0022001  Ba 

1.0018101  Ba 

1.0024242  Ba 

1.0010777  Ba 

1.0049429  Ba 

1.0024007  Ba 

1.0010719  Ba 

0.9999380  Ba 

1.0010719  B* 


-0.4258300 
-0.4281250 
-0.474442S 
-0.4281250 
-0.2941913 
-0.3997558 
-0.2480954 
-0.3487988 
-0.3038084 
-0.3740253 
-0.9492187 
-1.0715818 
-0.9038085 
-0.873583* 
-0.8721191 
-0.8971479 
-0.4405957 
-0.8805175 
-0.4084425 

2.2940912 
-0.4140444 
-0.4994428 
-0.2445800 
-0.3954054 

0.3308593 
-0.4279294 
-0.0544453 
-0.2947773 
-0.1212890 
-0.7454101 
-0.4801757 
-0.4208984 
-0.1928710 
-0.4208984 


THERMOCOUPLE   35 
THERMOCOUPLE   3( 
THERMOCOUPLE   37 
THERMOCOUPLE   38 
THERMOCOUPLE   39 
THERMOCOUPLE   40 
THERMOCOUPLE  41 
THERMOCOUPLE   42 
THERMOCOUPLE  43 
THERMOCOUPLE   44 
THERMOCOUPLE  45 
THERMOCOUPLE  44 
THERMOCOUPLE  47 
THERMOCOUPLE  48 
THERMOCOUPLE   49 
THERMOCOUPLE   50 
THERMOCOUPLE   51 
THERMOCOUPLE   52 
THERMOCOUPLE   53 
THERMOCOUPLE   54 
THERMOCOUPLE   55 
THERMOCOUPLE   54 
THERMOCOUPLE  57 
THERMOCOUPLE   58 
THERMOCOUPLE   59 
THERMOCOUPLE   60 
THERMOCOUPLE   61 
THERMOCOUPLE   62 
THERMOCOUPLE   63 
THERMOCOUPLE   64 
THERMOCOUPLE   45 
THERMOCOUPLE   64 
THERMOCOUPLE   67 
THERMOCOUPLE   68 


Ma 
Ma 


Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 


Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 
Ma 


1.0028639 

1.0017774 

1.0019427 

1.0019207 

1.0024245 

1.0024214 

1.0017195 

1.0013044 

1.0019274 

1.0024417 

0.9997830 

1.0019827 

1.0017757 

1.0030403 

1.0023384 

1.0024405 

1.0010144 

1.0013094 

1.0013409 

1.0013094 

1.0019S88 

1.0011234 

1.0017185 

1.0132902 

1.0026283 

1.0027475 

1.0014108 

1.0025404 

1.0019932 

1.0020733 

1.0018148 

1.0018148 

1.0018339 

1.0009155 


Ba 
Ba 
B> 
Ba 
Ba 
Ba 
Ba 
Ba 
Ba 
Ba 
Ba 
8* 
Ba 
Ba 
Ba 

Ba 
Ba 
Ba 
Ba 
Ba 
Ba 

Ba 

Ba 
Ba 

B* 

Ba 
Ba 
Ba 
Ba 
Ba 
Ba 
Ba 
Ba 
Ba 


-0.5545898 

-0.4930175 

-0.4770995 

-0.5316894 

-0.5033203 

-0.7974074 

-0.6256835 

-0.7514401 

-0.7330077 

-0.9266112 

-0.2326640 

-0.8003417 

-0.4931152 

-0.9011718 

-0.7045917 

-0.4196777 

-0.3534133 

-0.5520507 

-0.4940937 

-0.5520507 

-0.5768064 

-0.5678222 

-0.5253417 

-0.8309081 

-0.7043454 

-0.4194288 

-0.4912597 

-0.7730468 

-0.5812500 

-0.6718749 

-0.5980448 

-0.5980468 

-0.7208007 

-0.5204543 


Two  plots  of  T  actual  vs  Tmeasures  for  two  thermocouples  are 
presented  in  Figure  A.l  and  Figure  A. 2  to  demonstrate  their 
1  i  n  e  a  r  i  ty  . 
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APPENDIX  B 

TABULATED  DATA 

Heated  Cylinder  in  the  First  Row 


H  =  1.0 


Downstream  Row  No. 

1 

2 

3 

4 

5 
Heated    cylinder    in    the    first    row. 


1 . 5m/s 
Theta 
1.0 
.4  68 
.336 
.251 
.181 


H  =  4.6 


Downstream  Row  No. 

1 

2 

3 

4 

5 
Heated  cylinder  in  the  second  row. 


1 . 5m/s 
Theta 
1.0 
.497 
.365 
.276 
.197 


H  =  2.3 
Downstream  Row  No. 
1 
2 
3 
4 


1 . 5m/s 

Theta 

1.0 

.501 

.324 

.250 


3.0m/s 
Theta 
1.0 
.524 
.380 
.311 
.211 


3.0m/s 
Theta 
1.0 
.484 
.380 
.277 
.185 

3.0m/s 

Theta 

1.0 

.514 

.354 

.269 
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TABULATED  DATA 


Row 
No. 


s,o 


(°C) 


H  =  1.0 


H  =  2.3 


H  =  4.6 


1 
2 
3 
4 
6 

1 

2 
3 
4 
6 


1 
2 
3 
6 

1 
2 
3 
6 


1 
2 
3 
4 
6 

1 
2 
3 

4 
6 


98.5 
98.9 
99.6 
97.6 
96.3 

101.2 
89.4 
78.9 
96.8 
86.2 


111.0 

105.7 

102.8 

93.7 

99.3 
110.1 
118.2 

71.8 


96.7 

97.3 

104.8 


100 
91 

112 
99 
95 
91 


85.8 


(°C) 

Q(w) 

h(w/m2°c)      ReH 

Nu 

U    (m 

20.3 

3.03 

56.61 

394.2 

14.09 

1.5 

20.3 

3.05 

56.71 

394.2 

14.11 

1.5 

20.5 

3.28 

61.16 

393.9 

15.21 

i.5 

20.0 

3.18 

60.34 

394.5 

15.02 

1.5 

20.1 

3.74 

73.87 

394.4 

18.39 

1.5 

20.8 

4.59 

87.27 

787.2 

21.69 

3.0 

20.7 

4.17 

93.26 

787.4 

23.19 

3.0 

20.5 

3.73 

98.52 

787.9 

24.51 

3.0 

20.2 

4.74 

95.26 

787.6 

23.71 

3.0 

20.3 

4.57 

107.69 

788.3 

26.80 

3.0 

20.0 

4.12 

45.77 

907.4 

26.21 

1.5 

20.1 

4.04 

47.93 

907.1 

27.44 

1.5 

20.1 

4.13 

51.03 

907.1 

29.22 

1.5 

20.1 

3.78 

52.61 

907.1 

30.12 

1.5 

20.2 

5.50 

73.19 

1813.7 

41.90 

3.0 

20.1 

5.99 

69.86 

1814.2 

39.99 

3.0 

20.2 

6.81 

73.20 

1813.7 

41.90 

3.0 

20.1 

3.91 

80.03 

1814.2 

45.82 

3.0 

20.0 

3.27 

42.74 

1814.7 

48.96 

1.5 

20.2 

3.27 

42.49 

1813.7 

48.65 

1.5 

20.3 

3.71 

44.20 

1813.2 

50.60 

1.5 

20.6 

3.00 

36.9 

1811.6 

42.22 

1.5 

20.6 

2.92 

41.14 

1811.6 

47.06 

1.5 

19.5 

5.60 

62.83 

3  634.7 

72.04 

3.0 

19.8 

5.02 

65.70 

3631.6 

75.28 

3.0 

19.9 

4.80 

66.23 

3630.6 

75.87 
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